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PLASTIC DEFORMATIONS OF A BEAM UNDER IMPULSIVE LOADING! 


B. A. Cotter? and P. S. Symonds? 


ABSTRACT 


The deformation of a free-free beam subjected to a symmetric impulsive 
loading, specified by the initial velocity distribution of equation (1), is inves- 
tigated here by means of a plastic-rigid analysis. Bleich and Salvadori dis- 
cussed the same problem in a recent paper 1* 4 and presented results ob- 
tained from an elastic-plastic analysis and from a modified plastic-rigid 
analysis. In both analyses they assumed plastic hinge action to occur only 
at the center and that elsewhere the beam either deformed elastically or was 
rigid. The principal deformation magnitude is the final angle at the center 
of the beam. Results are presented here for this quantity in a general non- 
dimensional form, and comparisons are made with the corresponding results 
given by the methods used in 1*. The central angle determined by the com- 
plete plastic-rigid analysis differs markedly from that obtained by either of 
the “single-hinge” analyses. The range in which the complete plastic-rigid 
analysis can be expected to be valid is apparently different from that in which 
the single-hinge type of elastic-plastic analysis is accurate, the former giv- 
ing reasonably good results for the central angle only when the initial veloci- 
ties are sufficiently large, and the latter applying particularly to low veloci- 
ties. At low velocities the single-hinge type of plastic-rigid analysis and the 
complete plastic-rigid solution both involve appreciable errors (in opposite 
directions). 


INTRODUCTION 


The analysis given here is based on a method developed in a previous 
paper 2*. A moment-curvature relatior of the type shown in Fig. 1 is as- 
sumed to apply. Elastic deformations are neglected in comparison with 
plastic deformations. The curvature at any section is assumed constant 
until the moment magnitude reaches and maintains the limit moment Mp. 
Then the curvature may increase indefinitely, giving rise to a plastic hinge. 
Relative displacements of parts of the beam are assumed small enough so 
that changes in shape may be neglected. 

We consider a free-free uniform beam of length 24, limit moment M5, 
and mass per unit length m. The origin of the coordinate system is taken at 
the center of the beam with x denoting distance along the beam and y 


The results presented in this paper were obtained in the course of re- 
search sponsored by the Office of Naval Research under Contract 
N7onr-35801 with Brown University. 
Research Assistant, Division of Applied Mathematics, Brown University, 
Providence, R. I. 
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displacement, as shown in Fig. 2. The beam is given an initial symmetric 
velocity distribution as shown in Fig. 3 described by MT i¢ cos(F*) J 


where v, is a positive constant. The initial displacements are assumed zero. 
We are interested in finding the plastic deformation resulting from such an 
impulsive loading which is predicted by the plastic-rigid theory. In particu- 
lar the final angle of deformation at the middle of the beam is determined; 
this is an important measure of the damage to the beam. 


Plastic- Rigid Analysis 


After the initial impact the beam undergoes translation and deformation in 
four different stages before deformation is completed and the beam translates 
as a rigid body with constant, uniform velocity. Initially the beam may be 
described by 


y (x,o) = 0 ] 


(1) 
(x,0) = {1 + cos Ix] ¢ 2 if 


where y; stands for Oy/Ot. The boundary conditions for a free-free beam 
are that the shear force and bending moment vanish at the ends. Hence 


V(+ 2,t) 


M(+ 2,t) = t>0 


where V is the shear force and M is the bending moment. It can be seen from 
(1) that 


Yyx (%,0) = 0 <2 


Vor 


and (%,0) = (*,0) = - cos (ap) ial 22 


Figure 4 shows the initial rate of change of curvature along the beam. Hence 
initially the beam has zero curvature everywhere andfor O < Ix] < t 
the curvature is decreasing while for ‘ < |x| 4 4 the curvature is in- 
creasing. In view of the assumed moment-curvature relation of Fig. 1, and 
the fact that the moment and the rate of curvature must agree in sense, the 
initial moment distribution throughout the beam must be as shown in Fig. 5. 
In order to satisfy the boundary condition that M(+ dL, t) = O and to ensure 
that the bending moment is continuous along the beam, the moment distribu- 
tion at a somewhat later time t = t' must resemble that shown in Fig. 6. 


Phase 1 
Since the motion of the beam is symmetric about x = 0, we consider only 
the half beam defined by O < x < 4. Then at time t during the first phase 


| (2) 
i 
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the configuration of the half beam is that shown in Fig. 7a. In segment I the 
moment is constant and equal to -M,. Segments II and IV are rigid with 
angular velocities W; and Wo respectively measured positively in the clock- 
wise direction. Segment III has constant moment +M,. There are plastic 


hinges atx = nt, and where , and C are functions of time. 


E <0, n > O, and ¢ < O, where E stands for dE /dt, so that segments 
I and nr are diminishing in length. 
The velocity in segments I and III is simply 


y,(x,t) = cos(4*) } for x< 


since no shear forces act in these segments. In order for the transverse 
displacements to be continuous across the hinges, the transverse velocities 
must be continuous 2*. Hence 


y,(6%+4,t) = y, - « [1 + cos(nE ] 
y,(nt+,t) = (1 -,t) = [1 + cos(rn) 
y,Ct-,t) (1 + cost) 


It is evident that 


y,@t+,t) - (nt - 
| = 


Segment II, shown in Fig. 7b, undergoes no change in linear momentum dur- 
ing time t since it is not subjected to shearing forces from t = 0. Hence 
L(n-E) 


Cy, (ét+,t) -wyr Jar = 


Vo [ (ZX 
> 1 + } dx 
4 


Substituting the expression for y; (£4+,t) from (3) and that for w, from (4) 
in (5) and integrating, we obtain the identity 


(n-€)£ cos(n—) + cos(xn) J + [sin(xn) - sin(r—) = 
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External moments Mo act on the ends of this segment during the time inter- 
val zero tot. Hence an equation can be written expression the equality of 
angular impulse in the time t to the change of moment of momentum of II 


about the point 0. Then 


2M,t = Cy, -w rjr dr 


Jo 
[1 + cos (x2 )ax 


— the expressions for y, (E2+,t) and wy, as before and integrating we 
obtain 


M 
= (n-t)* + Jeos(nE) + [2(n-£)° - 5, (nn) 
6 


8 
(n-&) sin (xn). 


Now (6) may be rewritten in the following form with the aid of some well- 
known trigonometric identities: 


cos (BE +n) ] cos + sin = 0 
(9) 


From (9) we conclude that either 


cos[ B(E-n) + sin = 0 


cos [ B(E+n) (11) 


(10) yields & = n which corresponds to the solution at t = 0. From (11) the 
identity which must hold at t>0 is 


(12) 


Then (6) may be replaced by (12). 
Inserting (12) in (8) gives 


t= [(2n-1)° Jcos(xn) - S(2n-1) sin(nn) (13) 
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Equations (12) and (13) then give the location of two of the hinges during the 
first phase of the deformation. 

To determine the motion of the third hinge we observe from Fig. 7c that 
segment IV is acted on by an external moment M, at its end during the whole 
time t but is free of external forces. Since there is no change in linear mo- 
mentum 


L 
m Cy, t) worjdr = =| > {1 + 3X) ] dx (14) 
Jo 
cL 


Using the expression for y, (Ci+,t) in (3) and integrating we obtain an 
expression for Wo: 


{(1-C) cos (nC) + L (n%) (15) 
(1 


The equation for conservation of angular momentum, taking moments about 
P, becomes 


(ct, t) wp dr 


Va 


L 
-m {1 + cos (ZK) (x-C2L)dx (16) 


CL 
With the aid of (3), (16) reduces to 


M 
12 t= 3 + cos(n€) - + (14)3 + = 
° (17) 


Combining (15) and (17) gives 


= (1-¢)? - cos(nt) + (1-£) sin(xt) - 


MV (18) 

(18) then gives the location of the third hinge during the first phase. 
Expressions have thus been obtained for and as functions of 

and Ww, and Wp» as functions of —, n, and © during the first phase. It 

can be shown that & < 0, B > O and a < ©, so that sections I and III 

are diminishing in length. At sometime t,, (t,) = C(t,) and 

&(t,) > 0 so that III vanishes completely and the second phase of the mo- 


tion begins. If we define = E(t), 1) = V(t ,), and C, = 
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then fi» Ny, and 1 may be determined by observing that from (12) 
Setting t= ty, and = T in (13) and (17), and combining gives 


(2° = + sin(nn,) + =O (19) 


If (19) is solved numerically for ust the result is 


M 
&) = 0.2590 and ——2, ty = 2.69 x 10 
mv 


Phase 2 

For t > t, the motion of the beam corresponds to the configuration of 
Fig. 9a. Section I has constant moment -M,; sections II and IV are rigid 
with angular velocities w, andW,» respectively, positive when measured 
clockwise. Il and IV are separated by a moving hinge whose coordinate is 


x= nt where ) = (t). land IJ are joined at a moving interface at x = EL, 


where — = E(t), For this phase of the motion momentum equations can- 
not be written for II and IV separately since these sections have been sub- 
jected to unknown shearing forces and moments for t > ty. If we define 


21 = ¥p_(E44,t) and ao = yy,(ni+,t), then referring to Figs. 9b and 
9c, we can write the equations of Sotion for II and IV respectively as follows: 


m&n-£) Cay = (ne) 0 (20) 


ly = MB(1-9)? (21a) 


ay - - o (21b) 


aw 


2 
where and = Additional equations expressing the jump 


in transverse acceleration across the hinges at x = £2 and x = nt may be 
written, as discussed in 2*. 
The jump condition at the hinge at £2 is 


- = EL [O,(EL-,t) - 0,(EL4,¢)] 
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(EL+,t) = Wy and Yee (Et-, t) = 0 
Hence = sin(n—) - 
The jump condition at nf is 
But = EL - sin(n—) ] - 


and (nt -,t) =W) (nt+,t) Wo 


Hence 


Vo 


ag = sin(ng)] - Link) = (23) 


If (22) and (23) are inserted in (20b) and (21b), respectively, and the resulting 
expressions simplified the equations of motion become 


Yo 


(20, - sin(ng)) - 
(1-9)3 = 12 


(ME) + 2,405) + (1-n) = 0 


Thus we have to solve a set of four simultaneous first order non-linear dif- 
ferential equations for 7, & , w,, and w, as functions of t. A solution of 
this system can be obtained by means of numerical integration. First, how- 
ever, a simplification is made in which wy and wo are found in terms of 4, €, 
and t from momentum considerations and the resulting expressions substi- 
tuted in (24). 

The equation for conservation of linear momentum of the half beam for 
t > ty can be written 


M 

= - 24 
nd 
~)u, = 0 

(24) 
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Ain) 
Cy, t) - dr 


+ Cy, t) ~ wor 


v 
= {1 + ax 
EL 


If we observe that [1 + ) 
and y,(nt+,t) = -w 


then (25) simplifies to 
Vo(1-—) cos(m—) + = sin(n&) 


Since there is a moment -M, acting at x = 0 since t = 0 the equation for con- 
servation of angular momentum for the half beam, taking moments about 
x = 0, becomes 


A(n-€) 
"| - ar 


A 
+n (ni+r)dr 


+ cos(Z&)] x ax 
- m | = cos > x 
EL 


Completing the integration in (27) gives 


12 = 3(1-€2 2.) cos(n&) - 2 (1-7)? (24n) 
2 T Vo 


si 6 (28 
+ sin(xt) + % ) 


Equations (26) and (28) may be used to determine W, and wa in terms of £,7, 
and t and the results inserted in (24). If this is done and the following change 
of variables made 


M 
y=ln,2z2=14 and T = 6 


2 


(24) may be simplified to a pair of non-linear first order differential equa- 
tions in the dependent variables z and y only, namely 


dt ) z 


where (z) = 4 - 2°) cos(nz) + < z sin(nz) - + 


and = z cos(nz) + - 5) sin (nz) 


The differential equations (29) were integrated numerically using as initial 
conditions the values obtained for £, ", and t at the end of phase 1, namely 


2. = 1.614 x 1073 
2 
MV 


Initial values of y, and z, in the integration were obtained by the Kutta- 
Simpson rule applied to a system of two first-order differential equations 3*. 
Then the integration was completed over the required range by one of the 
standard methods of forward integration. The results of the integration show 
that 1 < 0 and € < O for this phase so that the interface atx = £4 and 
the hinge at x = nt are moving toward the center of the beam. At time 
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22° 

dt Cz-y le - (z-y) (29) 

yy = = 2590 = = 
=~ 


t = tg the angular velocities w and wz become equal while & is still positive, 
indicating the disappearance of the hinge atx = 1 4 . This concludes the 
second phase of the motion. If we define § (tg) = 9 and (tg) = Nog itis 
found that 


M 
= .1498 = and 2 = ,001216 
mv 


Phase 3 

For t > tg the beam consists of a rigid segment with angular velocity 
and a plastic segment with constant moment -M,. These are separated by a 
moving interface with coordinate x = £2 as shown in Fig. 10. A moment 
~Mbp still acts at x = 0. Momentum considerations may again be invoked to 
describe the motion of the beam in this phase. Since no external forces act 
on the half beam conservation of linear momentum gives 


p£(1-£) Lv 
m Cy, (et+,t) - wx) dx = ‘| {1 + cos(Z%) J dx 
EL 


(30) 


where y,(644,¢) = [1 + cos(n&) 


If the integration in (30) is carried out the result is 


Vo 


wh cos(nt) + 4 sin(né)] 


Conservation of angular momentum, taking moments about x = 0, gives 


(1-£ ) 
t "| - wx] (Er+x) dx 


1 + cos( ) < ax (32) 
With the value of w from (31), equation (32) reduces to 


M 
t = - (1-£)7] cos(n—) - +(1-£) sin(nE) + 

(33) 


(33) gives the location of the moving hinge for t > tg. It can be shown that 


E < O for this phase so that section I is diminishing in length. At some 
time tz, € = 0 and the half beam becomes entirely rigid with a hinge at x = 0. 
tz is simply determined by setting § = 0 in equation (33). Then 
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Mo — | 1 
t3 = - This concludes the third phase. 


For t > tg the half-beam is rigid and moves as if hinged at the point x = 0 
as shown in Fig. 11. Here the bending moment is -M) and if the angular ve- 
locity of the half-beam is w we may write the angular momentum equation, 
taking moments about x = 0, as follows: 


Cy,(o,t) - wx] x dx - =| {1 + cos(@X) ] x dx 


Mot = 
fe) (34) 


We observe that y;(o,t) is unknown for this phase. (34) reduces upon integrat- 
ing to 


M ¥+(0,t) 
MV Yo Vo ae (35) 


Conservation of linear momentum of the half-beam yields 
4 Vv 
=| Cy,(o,t) -wx)]dx =m > {1 + cos (BF) J dx (36) 


From (36) an expression is obtained for y;(o,t), namely 
Vo uw. 
y,(0,t) (1 + Ve (37) 


If (37) and (35) are combined the result is an expression for the angular ve- 
locity w as a function of time 


(38) 


The beam continues to translate and rotate with angular velocity w until 
some time t = tg when w = 0. From (38) 


(39) 


At t = tq deformation of the beam is completed and for t > tq the entire beam 
translates as a rigid body with constant velocity. The final transverse ve- 
locity vg of the beam results from setting w = 0 in (37) and is found to be = 


As a check, the same result is found directly by equating the initial linear 
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Phase 4 
i 
we -12 (4 t) 
Vo me my 


momentum mv, of the beam to the linear momentum in the final phase of 
rigid body translation, namely 2m) vy. 


Deformation at the Middle 


The angle of deformation at the middle of the beam may now be calculated 
quite readily. First observe that for 0<t <tg 


v 
¥,(x,t) = [1 + cos (5*)] x<¢ El 


v 
x 
andhence = y,,(x,t) = - 0,(x,t) = - > } 
Then the angle of deformation at the middle of the beam for 0 ¢ t < t 


is 
simply 6(0,t) = 0, For ts ty 3 


at t 
O(o4+,t) wdat' = (1 - 
L t, mv 


M M 
Then = - 6(—2, t- 


mv, mv 


Since at t = tg deformation ceases, the final angle 6(o+, ti) = Oof of de- 
formation at the middle is 


Mo 


Sor 


To summarize, then, if O(o+,t) = 9,(t), 


2 


t > ty 
M 
Figure 12 shows a plot of the non-dimensional quantity —— Oo(t) asa 
mv 


M ° 
fe) 
function of the non-dimensional time 


mv 


Energy Check 


An overall check on the results obtained in this analysis is available 
through energy considerations. The initial kinetic energy of the beam is 
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Mo t') 
|_| 
6.(t) = Sts ty (40) 


easily calculated and is found to be 3 miv,° . The final kinetic energy of 
the deformed beam is t miv°3 this is obtained from considerations of 


momentum conservation. Hence the loss of kinetic energy is g miv,°. 


This must be equal to the plastic work done on the beam during deformation, 
Wo» which is given by 


= - 


where O,¢ = and = O(L, te). One and may be 


found by direct integration over the entire range O < t <¢ t» with the ex- 
ception of phase 2 during which the integration must be performed numerical- 


ly. If this is done the result is On¢ = .0621 Wo and ©. = 
2 
Then Wy = .1252 which agrees within the 


accuracy of the calculations with the loss in kinetic energy. 


Discussion of Results 


We observe that, since all deformation at the midpoint takes place during 
phase 4, @_ (t) may be calculated from the motion of this phase once the 
value of tg has been obtained from the solution of phase 3. Hence it appears 
that 6 of may be determined without obtaining a solution for the first two 
configurations thereby obviating the need for resorting to numerical integra- 
tion. However, in order to ensure that no phases of the type shown in Fig. 

13a and 13b occur in place of 2 and 3, it is necessary to determine the motion 
for these two phases. If the configurations of Fig. 13 were to apply the cen- 
tral angle @ ,.(t) would not be easily calculable. 

A comparison may be made between the present results and those obtained 
by Bleich and Salvadori in 1*. They obtain midpoint angles of deformation for 
a specific beam with m = .800, 4 = 160, My = 6.69 x 104, and EI = 3.58 x 108 
for two values of vo, namely v, = 11.1, and v, = 67.4. Their results for the 
two cases appear in Table I together with values of 6 9(t) obtained by the cor- 
rect rigid-plastic analysis developed here. 6, is obtained in 1* by an elastic- 
plastic analysis in which the limit moment Mb) is assumed to occur only at 
the mid-point at an “initial plasticity time,” to. 9, is the result determined 
in 1* by assuming the beam hinged at the middle and rigid elsewhere from 
t = 0 until deformation is completed. 6) is obtained from the same “single- 
hinge” assumption for t >to with the beam elastic for 0 <t< ty. Since the 
angles of deformation reported in 1* are defined to be equal to the discontinu- 
ity in slope: 


= y,(o+,t) -y,.0-,t) 


the quantities 6 ,, Oy , and Ge have, by definition, twice the magnitude of the 
central angles discussed in the present analysis. Hence in Table I we com- 
pare 265 with Oy, 6), and 6¢. 
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If we employ the criterion for the validity of the plastic-rigid analysis 
suggested in 2*, namely that the plastic work done on the beam, be be much 
25 
greater than the maximum elastic energy te, then, if R = —e— our 
I 
results should be validforR >> 1. For vo = 11.1, R = .986 and for vo = 
67.4, R = 36.3. Thus we might expect 6 of to be approximately correct for the 
higher value of v, only. However, Table I indicates marked differences be- 
tween @, and the final angles obtained by either of the “single-hinge” analy- 
ses for Vg = 67.4; the final angle of deformation which is predicted by the 
three different analyses in 1* is about 50% higher than that yielded by the 
rigid-plastic solution. For vg = 11.1, R ~ 1 so that the value of 6 of obtained 
in this case is not expected to be valid. However, the results of the single- 
hinge elastic-plastic solution are probably accurate here. For this velocity 
the positive error of the single-hinge plastic-rigid solution has about the 
same magnitude as the negative error of the complete plastic-rigid solution. 

It should be noted that the disagreement between the elastic-plastic and 
the present rigid-plastic solution at the higher velocity vo - 67.4 does not 
mean that the latter solution must be in error. In fact, the agreement be- 
tween the elastic-plastic solution and the single-hinge rigid-plastic solution 
is an indication that the elastic-plastic solution is inaccurate at this velocity 
because it takes account of plastic flow solely at the mid-section. Such 
agreement between the two “single-hinge” solutions is to be expected at suf- 
ficiently high velocities, if these are regarded as the correct solution for a 
problem of a fictitious beam whose yield stress varies along the beam so that 
plastic flow at the limit moment Mb occurs only at the middle. Both solutions 
are presumably incorrect for the uniform beam problem. The complete 
plastic-rigid analysis of this paper may well be accurate at the higher veloci- 
ty, but confirmation by comparison with a complete elastic-plastic solution 
still remains to be obtained. 

In an oral communication based on an advance copy of the present paper 
Professor Bleich has called attention to the interest of a comparison between 
the final angle at the mid-point determined by the elasto-plastic analysis of 
1* and the total angle across the central plastic zone in the complete plastic- 
rigid analysis. In the latter solution there is a zone of length { with negative 
permanent curvature occupying the middle portion of the beam, while the two 
outer quarters have positive curvature. The total angle between the tangents 


at the two ends of the negative zone is 26 0.1242 » rer 
the two cases treated numerically in 1* the final values of the relative angle 
across the entire negative plastic zone of length! are: 


Vo = 


29 
x 10° = 1.60. 
Vo 


The very close agreement of these results with the corresponding values 
(2.63, 1.58) of the single-hinge rigid-plastic solution is natural in view of 

the fact that the work done in positive plastic deformation (in the outer quar- 
ters of the beam) is less than one percent of the work done in producing nega- 
tive plastic deformation; this is shown by the figures given on p. 17. 
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at 
Vo = lel 
20 
Vo 


As far as damage to the beam is concerned the angle of slope discontinuity 
at a plastic hinge is especially significant. At any plastic hinge there must of 
course be a small but finite plastic zone in a physical beam, since a slope 
discontinuity would imply infinite strains. Such a smoothing out of the theo- 
retical angle of discontinuity would occur at the central hinge predicted either 
by the elasto-plastic analysis of 1* or by the present plastic-rigid analysis. 
The comparisons between central angles shown in Table I appear therefore 
to be of primary interest in assessing the various methods. 
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Table I 


ANGLE OF DEFORMATION AT MID-POINT OF FREE-FREE BEAM 
Vo = 11.1 Vo = 67.4 


Vo Vo Vo Vo 
0 0 0 


- Correct Plastic-Rigid Solution 


- Upper Bound Solution 


- Elastic-Plastic Solution | Bleich and 
Salvadori 

- Lower Bound Solution 

~ 


| 
| t 200x103 || Se 02) 
Yo Vo Vo Vo 
0 0 | 
205 66 9) 9) | 025 
-061 0 73 | 0 
.10 | .450 | 1.30 0 421 93 | 
-932 | 1.80 0 1.0 1.16 | 
20 1.310 | 2.18 0 1.25 818 1.32 | 
25 | 1.574% | 2.44 111.50 1.46 
30 1.72% | 2.57 | 1.38 .68 {11.75 | 1.035 1.53 
35 1.770 | 2.63 | 1.97 | .8% |/2.0 | 1.073] 1.56 
I-40 2.63 | 2.20 | .93 [l2.25! 1.075 | 1.58 
| | £95 | | 
| 

| 
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PROCEEDINGS-SEPARATES 


The technical papers published in the past year are presented below. Technical-division 
sponsorship is indicated by an abbreviation at the end of each Separate Number, the symbols 
referring to: Air Transport (AT), City Planning (CP), Construction (CO), Engineering Mechanics 
(EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), Sanitary 
Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping 
(SU), and Waterways (WW) divisions. For titles and order coupons, refer to the appropriate 
issue of “Civil Engineering” or write for a cumulative price list. 


VOLUME 80 (1954) 
APRIL: 428(HY)°, 429(EM)°, 430(ST), 431(HY), 432(HY), 433(HY), 434(ST). 
MAY: 435(SM), 436(CP)°, 437(HY)°, 438(HY), 439(HY), 440(ST), 441(ST), 442(SA), 443(SA). 


JUNE: 444(SM)°, 445(SM)®, 446(ST)©, 447(ST)®, 448(ST)®, 449(ST)®, 450(ST)®, 451(ST)®, 452(SA)®, 
453(SA)®, 454(SA)®, 455(SA)®, 456(SM)¢. 


JULY: 457(AT), 458(AT), 45S(AT)°, 460(IR), 461(IR), 462(IR), 463(1R)°, 464(PO), 465(PO)°. 


AUGUST: 466(HY), 467(HY), 468(ST), 469(ST), 470(ST), 471(SA), 472(SA), 473(SA), 474(SA), 
475(SM), 476(SM), 477(SM), 478(SM)°, 479(HY)°, 480(ST)©, 481(SA)°, 482(HY), -483(HY). 


SEPTEMBER: 484(ST), 485(ST), 486(ST), 487(CP)°, 488(ST)°, 489(HY), 490(HY), 491(HY)°, 
492(SA), 493(SA), 494(SA), 495(SA), 496(SA), 497(SA), 498(SA), 499(HW), 500(HW), 501(HW)°, 
502(WwW), 503(WW), 504(Ww)°, 505(CO), 506(CO)®, 507(CP), 508(CP), 509(CP), 510(CP), 
511(CP). 


OCTOBER: 512(SM), 513(SM), 514(SM), 515(SM), 516(SM), 517(PO), 518(SM)°, 519(IR), 520(IR), 
521(IR), 522(IR)°, 523(AT)©, 524(SU), 525(SU)©, 526(EM), 527(EM), 528(EM), 529(EM), 
530(EM)°, 531(EM), 532(EM), 533(PO). 


NOVEMBER: 534(HY),535(HY), 536(HY), 537(HY), 538(HY)°, 539(ST), 540(ST), 541(ST), 542(ST), 
543(ST), 544(ST), 545(SA), 546(SA), 547(SA), 548(SM), 549(SM), 550(SM), 551(SM), 552(SA), 
553(SM)°, 554(SA), 555(SA), 556(SA), 557(SA). 


DECEMBER: 558(ST), 559(ST), 560(ST), 561(ST), 562(ST), 563(ST)°, 564(HY), 565(HY), 566(HY), 
567(HY), 568(HY)°, 569(SM), 570(SM), 571(SM), 572(SM)°, 573(SM)°, 574(SU), 575(SU), 576(SU), 
577(SU), 578(HY), 579(ST), 580(SU), 581(SU), 582(Index). 


VOLUME 81 (1955) 


JANUARY: 583(ST), 584(ST), 585(ST), 586(ST), 587(ST), 588(ST), 589(ST)©, 590(SA), 591(SA), 
592(SA), 593(SA), 594(SA), 595(SA)°, 596(HW), 597(HW), 598(HW)°,599(CP), 600(CP), 601(CP), 
602(CP), 603(CP), 604(EM), 605(EM), 606(EM)°, 607(EM). 


FEBRUARY: 608(WW), 609(WW), 610(WW), 611(WW), 612(WW), 613(WW), 614(WW), 615(WW), 
616(WW), 617(IR), 618(IR), 619(01R), 620(IR), 621(IR)°, 622(IR), 623(IR), 624(HY)°, 625(HY), 
626(HY), 627(HY), 628(HY), 629(HY), 630(HY), 631(HY), 632(CO), 633(CO). 


MARCH: 634(PO), 635(PO), 636(PO), 637(PO), 638(PO), 639(PO), 640(PO), 641(PO)*, 642(SA), 
643(SA), 644(SA), 645(SA), 646(SA), 647(SA)°, 648(ST), 649(ST), 650(ST), 651(ST), 652(ST), 
653(ST), 654(ST)©, 655(SA), 656(SM)°, 657(SM)°, 658(SM)°. 


APRIL: 659(ST), 660(ST), 661(ST)°, 662(ST), 663(ST), 664(ST)©, 665(HY)°, 666(HY), 667(HY), 
668(HY), 669(HY), 670(EM), 671(EM), 672(EM), 673(EM), 674(EM), 675(EM), 676(EM), 677(EM), 
678(HY). 


c. Discussion of several papers, grouped by Divisions. 
e. Presented at the Atlantic City (N.J.) Convention in June, 1954. 
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